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Abstract 

We study the soliton type solutions arising in two-dimensional quantum 
chromodynamics (QCD2). In bosonized QCD2 these type of solutions emerge 
as describing baryons and quark solitons (excitations with "colored" states), 
respectively. The so-called generalized sine-Gordon model (GSG) arises as the 
low-energy effective action of bosonized QCD2 for unequal quark mass parame- 
ters, and it has been shown that the relevant solitons describe the normal and 
exotic baryonic spectrum of QCD 2 [JHEP(03)(2007)(055)]. In the first part of 
this chapter we classify the soliton and kink type solutions of the sl(3) GSG 
model with three real fields, which corresponds to QCD2 with three flavors. 
Related to the GSG model we consider the sl(3) affine Toda model coupled to 
matter fields (Dirac spinors) (ATM). The strong coupling sector is described 
by the si (3) GSG model which completely decouples from the Dirac spinors. In 
the spinor sector we are left with Dirac fields coupled to GSG fields. Based on 
the equivalence between the U(l) vector and topological currents, which holds 
in the ATM model, it has been shown the confinement of the spinors inside the 
solitons and kinks of the GSG model providing an extended hadron model for 
"quark" confinement [JHEP(01)(2007)(027)]. Moreover, it has been proposed 
that the constituent quark in QCD is a topological soliton. These qualitons 
(quark solitons), topological excitations with the quantum numbers of quarks, 
may provide an accurate description of what is meant by constituent quarks 
in QCD. In the second part of this chapter we discuss the appearance of these 
type of quark solitons in the context of bosonized QCD 2 (with Nf — 1 and 
N c colors) and the relevance of the si (2) ATM model in order to describe the 
confinement of the color degrees of freedom. We have shown that QCD 2 has 
quark soliton solutions if the quark mass is sufficiently large. 



1 Introduction 

A useful theoretical laboratory for studying several problems in Quantum Chromo- 
dynamics is QCD in two dimensions [1, 2]. This theory can be written in bosonized 
form [3] for arbitrary numbers of colors N c and flavors Nf [4]. It reflects accurately 



1 



Baryons and quark solitons in QCD 2 



2 



the phenomena of quark confinement and condensation in the vacuum that we ex- 
pect to occur in QCD in four dimensions. In the low-energy and strong coupling 
limit (e c >> m q , e c =coupling constant, m g =quark mass ) QCD2 has finite-energy 
soliton solutions for arbitrary values of N c and Nf that can be interpreted as baryons 
[2], in close analogy with the skyrmion interpretation of baryons as solitons in QCD4 
[5]. In this limit the static classical soliton which describes a baryon in QCD2 turns 
out to be the ordinary sine-Gordon (SG) soliton. It has been shown that various 
aspects of the low-energy effective QCD2 action with unequal quark masses can be 
described by the so-called (generalized) sine-Gordon model (GSG) [6]. 

Moreover, it has been proposed that the constituent quark in QCD4 is a topo- 
logical soliton [7]. These qualitons (quark solitons), topological excitations with 
the quantum numbers of quarks, may provide an accurate description of what is 
meant by constituent quarks in QCD. Related to this phenomenon, it has been 
found certain static soliton solutions to QCD2 that have the quantum numbers of 
quarks [8]. They exist only for quarks heavier than the dimensional gauge coupling 
(e c << m q ), and have infinite energy, corresponding to the presence of a string 
carrying the non-singlet color flux off to spatial infinity. 

On the other hand, the sine-Gordon model (SG) has been studied over the 
decades due to its many properties and mathematical structures such as integrability 
and soliton solutions. It can be used as a toy model for non-perturbative quantum 
field theory phenomena. In this context, some extensions and modifications of the 
SG model deserve attention. An extension taking multi-frequency terms as the 
potential has been investigated in connection to various physical applications [9, 
10, 11, 12]. Another extension defined for multi-fields is the so-called generalized 
sine-Gordon model (GSG) which has been found in the study of the strong/weak 
coupling sectors of the so-called sl(N,<E) affine Toda model coupled to matter fields 
(ATM) [14, 15, 16]. In connection to these developments, the bosonization process 
of the multi-flavor massive Thirring model (GMT) provides the quantum version 
of the (GSG) model [17]. The GSG model provides a framework to obtain (multi- 
) soliton solutions for unequal mass parameters of the fermions in the GMT sector 
and study the spectrum and their interactions. The extension of this picture to the 
NC space-time has been addressed (see [18] and references therein). 

It has been conjectured that the low-energy action of QCD2 (e >> m q , m q quark 
mass and e gauge coupling) might be related to massive two dimensional integrablc 
models, thus leading to the exact solution of the strong coupled QCD2 [2]. In 
particular, it has been shown that the sl(2) ATM model describes the low-energy 
spectrum of QCD2 (1 flavor and N c colors) and the exact computation of the string 
tension was performed [19]. A key role has been played by the equivalence between 
the Noether and topological currents at the quantum level. Moreover, one notice 
that the SU(n) ATM theory [14, 15] is a 2D analogue of the chiral quark soliton 
model proposed to describe solitons in QCD4 [20] , provided that the pseudo-scalars 
lie in the Abelian subalgebra and certain kinetic terms are supplied for them. 
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Besides, coupled systems of scalar fields have been investigated by many au- 
thors [21, 22, 23, 24, 25, 26]. One of the motivations was the study of topological 
defects in relativistic field theories; since realistic theories involve more than one 
scalar field, the multi-field sine-Gordon theories with kink-type exact solutions de- 
serve some attention. The interest in the study of the classical limit of string theory 
on determined backgrounds has recently been greatly stimulated in connection to 
integrability. It has been established that the classical string on R x S 2 is essentially 
equivalent to the sine-Gordon integrable system [27]. More recently, on Rx S 3 back- 
ground utilizing the Pohlmeyers reduction it has been obtained a family of classical 
string solutions called dyonic giant magnons which were associated with solitons 
of complex sine-Gordon equations [28]. String theory on R x S N ~ 1 is classically 
equivalent to the so-called SO(N) symmetric space sine-Gordon model (SSG) [29]. 

In the first part of this chapter we study the spectrum of solitons and kinks 
of the GSG model proposed in [14, 15, 17] and consider the closely related ATM 
model from which one gets the classical GSG model (cGSG) through a gauge fixing 
procedure. Some reductions of the GSG model to one-field theory lead to the usual 
SG model and to the so-called multi-frequency sine-Gordon models. In particular, 
the double (two-frequency) sine-Gordon model (DSG) appears in a reduction of the 
sZ(3,C) GSG model. The DSG theory is a nonintegrable quantum field theory with 
many physical applications [11, 12]. 

In the ATM model, once a convenient gauge fixing is performed by setting to 
constants some spinor bilinears, we are left with two sectors: the cGSG model which 
completely decouples from the spinors and a system of Dirac spinors coupled to the 
cGSG fields [16]. In the references [30, 31] a 1 + 1-dimensional bag model for quark 
confinement is considered, we follow their ideas and generalize for multi-flavor Dirac 
spinors coupled to cGSG solitons and kinks. The first reference considers a model 
similar to the sl(2) ATM theory, and in the second one the DSG kink is proposed 
as an extended hadron model. 

In the second part of this chapter we examine the quark soliton type solutions 
in QCD2. Regarding this phenomenon several properties of the ATM model deserve 
careful consideration in view of the relationships with two-dimensional QCD. For 
simplicity we concentrate on the si (2) ATM model. So, in order to disentangle the 
quark solitons one needs to restore the heavy fields, i.e. the fields associated to 
the color degrees of freedom. This is done in two steps. First, by including N c 
dynamical Dirac spinors coupled to the Toda field, second by breaking the chiral 
symmetry through certain bilinear terms in the scalar fields of the bosonized effective 
Lagrangian. In this way we arrive at a model similar to the one proposed in [8] in 
the regime when m q » e c . We have shown that QCD2 has quark soliton solutions 
if the quark mass is sufficiently large. 

In the next section we define the sZ(3) GSG model and study its properties such 
as the vacuum structure and the soliton, kink and bounce type solutions. In section 
3 we consider the sZ(3) affine Toda model coupled to matter and obtain the cGSG 
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model through a gauge fixing procedure. It is discussed the physical soliton spectrum 
of the gauge fixed model. In section 4 the topological charges are introduced, as well 
as the idea of baryons as solitons (or kinks) , and the quark confinement mechanism 
is discussed. In section 5 we examine the quark soliton solutions of QCD2 and 
discuss the role played by the effective sl(2) ATM model. The discussion section 
outlines the main results of this contribution and some lines of future research. In 
appendix A we provide the zero curvature formulation of the si (3) ATM model. 



2 The GSG model 

The generalized sine-Gordon model (GSG) related to sl(N) is defined by [14, 15, 17] 

Nf 

S= [d 2 xJ2 [^(^) 2 + ^(cosA^-l)]. (1) 

i=l 

The <&i fields in (1) satisfy the constraints 

$ p =J2<r P i*i, p = N,N + l,...,N f , N f = N{N ~ 1 \ (2) 
i=i 

where a p i are some constant parameters and Nf is the number of positive roots of 
the Lie algebra sl(N). In the context of the Lie algebraic construction of the GSG 
system these constraints arise from the relationship between the positive and simple 
roots of sl(N). Thus, in (1) we have (N — 1) independent fields. 

We will consider the sl(3) case with two independent real fields 931,2, such that 

$1 = 2^1-^2; $2 = 2(^2 - <Pi ; $3 = r ipi + s (p 2 , s,rGlR (3) 
which must satisfy the constraint 

/?3$ 3 = Sxfr&l + 5 2 P2®2, A = AM, (4) 

where /3b, fj, 61,82 are some real numbers. Therefore, the sl(3) GSG model can be 
regarded as three usual sine-Gordon models coupled through the linear constraint 
(4). 

Taking into account (3)-(4) and the fact that the fields cp± and (f2 are independent 
we may get the relationships 

I 2s + r 

V282 = ^3 = — ; — (Mi + V282); po = - — ■ — (5) 

r + s 2r + s 

The sl(3) model has a potential density 

3 

^1 = ^^(1-008/3^) (6) 
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Figure 1: GSG potential V for the parameter values ui = 1/2, 5\ = 2, 5 2 = 
1, f 2 = 1, r = s = 1, A) = 1, Ml = y«2 = 1- 



The GSG model has been found in the process of bosonization of the generalized 
massive Thirring model (GMT) [17]. The GMT model is a multiflavor extension 
of the usual massive Thirring model incorporating massive fermions with current- 
current interactions between them. In the sl(3) construction of [17] the parameters 
Si depend on the couplings /3, and they satisfy certain relationship. This is obtained 
by assuming [i{ > and the zero of the potential given for $j = -£rii, which 
substituted into (4) provides 

niSi + n 2 5 2 = n 3 , rij G 2 (7) 

The last relation combined with (5) gives 

( 2r + s )^ + (2 S + r)^ = 3^. (8) 
v\ u 2 1/3 

The periodicity of the potential implies an infinitely degenerate ground state and 
then the theory supports topologically charged excitations. A typical potential is 
plotted in Fig. 1. The vacuum configuration is related to the fundamental weights 
(see sections 3, 4 and the Appendix). For the moment, consider the fields $1 and 
$2 and the vacuum lattice defined by 

(* l ,* 2 ) = ^,^) ) n a EZ. (9) 
PO V\ v 2 

It is convenient to write the equations of motion in terms of the independent 
fields ip\ and ip 2 

d 2 ip\ = -/zi/5iAnsin[/3i(2<pi - (p 2 )] - fj, 2 f3 2 Ai 2 sm[f3 2 (2(p 2 - <pi)] + 

// 3 /?3Ai3sin[/3 3 (7-931 + stp 2 )] (10) 

d 2 ip 2 = -/Ji/3iA 2 isin[/3i(2v9i - ip 2 )] - fi 2 {3 2 A 22 sm[f3 2 (2ip 2 -<pi)) + 

M3/33A 2 3sin[/3 3 (rv9i + stp 2 )], (11) 
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where 



A = p 2 u 2 (4 + 5 2 + 6 2 p 2 r 2 ), B = f3 2 v 2 {l + A5 2 + S 2 p 2 s 2 ), 

C = $v 2 (2 + 25 2 + S 2 lP 2 1 rs), 

An = (C-2B)/A, A 12 = (B - 2C)/A, A 13 = (r B + sC)/A, 

A 21 = (A-2C)/A, A 22 = (C - 2A)/A, A 23 = (r C + s A) /A 




Notice that the eqs. of motion (10)-(11) exhibit the symmetry 



(12) 



Some type of coupled sine-Gordon models have been considered in connection to 
various interesting physical problems [32]. For example a system of two coupled SG 
models has been proposed in order to describe the dynamics of soliton excitations in 
deoxyribonucleic acid (DNA) double helices [33]. In general these type of equations 
have been solved by perturbation methods around decoupled sine-Gordon exact 
solitons. 

The system of equations (10)-(11) for certain choice of the parameters r and s 
will be derived in section 3 in the context of the sl(3) ATM type models, in which the 
fields ifi and <£> 2 couple to some Dirac spinors in such a way that the model exhibits 
a local gauge invariance. The ATM relevant equations of motion have been solved 
using a hybrid of the Hirota and Dressing methods [34]. However, in this reference 
the physical spectrum of solitons and kinks of the theory, related to a convenient 
gauge fixing of the model, have not been discussed, even though the topological and 
Noether currents equivalence has been verified. The appearance of the so-called tau 
functions, in order to find soliton solutions in integrable models, is quite a general 
result in the both Dressing and Hirota approaches. In this section, we will find 
soliton and kink type solutions of the GSG model (10)-(11) and closely follow the 
spirit of the above hybrid method approach to find soliton solutions. 

The general tau function for an n— soliton solution of the gauge unfixed ATM 
model has the form [34, 35] 



2 



r 





(13) 



Pl Pn=0 



Zi = ji(x-Vit), c pi p „eC 



Since the GSG model describes the strong coupling sector (soliton spectrum) 
of the ATM model [14, 15] then one can guess the following Ansatz for the tau 
functions of the GSG model 
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where the tau functions Tj (i = 0,1,2) are assumed to be of the form (13). We 
will see that the Ansatz (14) provides soliton and kink type solutions of the model 
(10)-(11), in this way justifying a posteriori the assumption made for the relevant 
tau functions. 

Assuming that the fields ip a (a = 1, 2) are real, from (14) one can write 



Vi, 2 



arctan[F(ri i2 ,r ) 



(15) 



F = 



ei 



ei([i?e(r lj2 )] 2 + [/m(r li2 )] 2 ) - (&(n, 2 )Re(T ) + Im{r lt 2 )/m(ro)) 



7m(ri j2 ) * i?e(ro) - i2e(ri )2 ) * Im(T ) 



±1 



(16) 



In terms of the tau functions the system of equations (10)-(11) becomes 



2t 



0§ 



^ 2 n {dnf d 2 r , (dr ) 2 i , /? im A n r(r 2 r ) 4 ^ -r^ 1 



Tl 



TO 



+ 



+ 



2i 



(t 2 t ) 



2j/i t-4i/i 



+ 



A 12 r (rin,) 4 ^ - rf* 1 /3 3 ^ 3 A 13 r (t ) 4 ^h-*) - r^r^ 3 



2i I (nro) 2 ^^ 
3 2 r 2 (<9r 2 ) 2 a 2 r 



2t ^ 2 
r 2 

/3 2 /i 2 A 22 
2z 



(5r ) 



2i L (r2) 2^a( ri )2r, 3T 2 ^(^) 

toA 21 r(r 2 r ) 4 ^ - 1 



0, (17) 



to 



+ 



(TlT 



1 4^2 



L (Tir ) 2 ^T 2 4l/2 



2* L (r.ro) 2 ^^ 1 

/33^3A23r(To) 4 ^^ +S )-T 1 4 ^r 2 4 ^ 



+ 



2i L (r2) 2s, 3 ( ri )2r,3 T 2 ^(^) 



0. (18) 



We will see that the 1-soliton and 1-kink type solutions are related to half- 
integer or integer values of the parameters i/j and the values r,s = 0,1. In the next 
subsections we write the 1-antisoliton, 1-antikink and bounce type solutions, and in 
order to perform the cumbersome computations we resort to the MAPLE program. 

2.1 One soliton associated to <pi 

Consider the tau functions 

ro = 1 + i dexp[y(x — vt)]; t\ = 1 — i dexp[j(x — vt)]; r 2 = 1 + idexp[j(x — vt)]. 
This choice satisfies the system of equations (17)-(18) for the set of parameters 
v x = 1/2, «5i = 2, S 2 = 1, v 2 = 1, u z = 1, r = 1. (19) 

provided that 

1 



13/x 3 = 5/i 2 - 4/xi, 7 X = ^(6yU2 + 3/ii). 



(20) 
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Now, taking e\ = 1 in Eq. (16) and the relation (15) one has 
4 

tpi = — — arctan{<i exp[7i(x — vt)]}, ip 2 = 0. (21) 
Po 

This solution is precisely the sine-Gordon 1-antisoliton associated to the field ipi 
with mass M\ = We plot a soliton of this type in Fig. 3. 

2.2 One soliton associated to p 2 

Next, let us consider the tau functions 

ro = 1 + i dexp[j(x — vt)], t\ = 1 + idexp[j(x — vt)], t 2 = 1 — i dexp[y(x — vt)] 
This set of tau functions solves the system (17)-(18) for the choice of parameters 
v x = 1, <5i = 1, 5 2 = 2, v 2 = 1/2, 1/3 = 1, 8 = 1 (22) 

provided that 

13^3 = 5^1-4^2, 72 = 7§( 6/Ul + 3 ^) (23) 

Now, choose e\ = 1 in (16) and through (15) one can get 
4 

ip 2 = — — arctan{dexp[7 2 (x - vt)]}, = (24) 
Po 

Similarly, this is the sine-Gordon 1-antisoliton associated to the field (p 2 with 
mass M 2 = ^ and its profile is of the type shown in Fig 3. 

2.3 Two one-solitons associated to (p = ^ lj2 
Now, let us consider the tau functions 

To = 1 + i d exp[7(x — vt)], t± = 1 — i d exp[j(x — vt)], t 2 = 1 — i d exp[y(x — vt)]. 
This choice satisfies (17)-(18) for 

Ul = 1, S 1 = 1/2, v 2 = 1, 6 2 = 1/2, v 3 = 1/2, r = s = l, (25) 
provided that 

d 2 = 1, 387! = 25^1 + 13^2 + 19^3 (26) 
Now, taking e\ = 1 in (15) one has 

<Pi = <P2 = <Pi, (27) 
4 

<£i = — — arctanjd exp[73(x — vt)]}. (28) 
Po 
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This is a sine-Gordon 1-antisoliton associated to both fields ip\^ in the particular 

case when they are equal to each other. It possesses a mass M3 = ^S-. 

Po 

In view of the symmetry (12) we are able to write 

d 2 = l, 3871 = 25^ + 13^1 + 19^3, (29) 
and then on has another soliton of this type 

Pi = ^2 = 02, (30) 

4 

02 = — — arctanjd exp[74(x — vt)]}. (31) 
Po 

It possesses a mass M4 = ^M-. This 1-antisoliton is of the type shown in Fig. 3. 

The GSG system (10)-(11) reduces to the usual SG equation for each choice 
of the parameters (19), (22) and (25), respectively. Then, the n— soliton solutions 
in each case can be constructed as in the ordinary sine-Gordon model by taking 
appropriate tau functions in (13)-(14). 

The baryon number associated to each of the above 1-soliton solutions has been 
computed in connection to QCD2, and it takes the same value B = N c (in this 
normalization the quark has baryon number B quar k = 1) [6]. 

A modified model with rich soliton dynamics is the so-called stepwise sine- 
Gordon model in which the system parameter depends on the sign of the SG field 
[36]. It would be interesting to consider the above GSG model along the lines of 
this reference. 

2.4 Mass splitting of solitons 

It is interesting to write some relations among the various soliton masses 

Mf = -^(109Mf + 5Mf); M 4 2 = ^(109M? + 5M|); (32) 

If /Hi = fJ-2 then we have the degeneracy M\ = M 2 , and M3 = M4 = -y/3/ 2M\ . 
Notice that if Mi / M 2 then M 3 < Mi + M 2 and M 4 < M x + M 2 , and the third and 
fourth solitons are stable in the sense that energy is required to dissociate them. 

2.5 Kinks of the reduced two-frequency sine-Gordon model 

In the system (10)-(11) we perform the following reduction (p = ipi = (p 2 such that 

$1 = $2, $3 = 9$1, (33) 

with q being a real number. Therefore, using the constraint (4) one can deduce the 
relationships 
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Moreover, for consistency of the system of equations (10)-(11) we have to impose 
the relationships 

I/l/XlAn + ^ 2 /"2Ai2 = I>lAtlA 2 l + ^2/"2A 2 2, (35) 

Ai 3 = A 23 . (36) 

These relations imply 

5 2 = 1, fn = 5fi 2 (37) 



Taking into account the relations (34) and (37) together with (5) we get 

Mi =1*2, S = 1, y\ = v 2 , ^3 = r = s = 1. (38) 

Thus the system of Eqs.(10)-(11) reduce to 

<9 2 <I> = _^l sin ( I/l$ )_ ^igin^^^), $ = (39) 

This is the so-called two-frequency sine-Gordon model (DSG) and it has been 
the subject of much interest in the last decades, from the mathematical and physical 
points of view. It encounters many interesting physical applications, see e.g. [11, 
12, 31, 32]. 

If the parameter q satisfies 

71 

? = -eQ (40) 

m 

with m, n being two relative prime positive integers, then the potential ^%(1 — 



cos(^i<I>)) + ^(1 — cos{qv\<&)) associated to the model (39) is periodic with period 

— m = n. (41) 

As mentioned above the theory (39) possesses topological excitations. The fun- 
damental topological excitations degenerates in the >x\ = limit to an n— soliton 
state of the relevant sine-Gordon model and similarly in the limit ^3 = it will be 
an m-soliton state. For general values of the parameters fi±, n.3, Si, v\ the solitons 
are in some sense "confined" inside the topological excitations which become in this 
form some composite objects. On the other hand, if q ^ Q then the potential is 
not periodic, so, there are no topologically charged excitations and the solitons are 
completely confined [9, 10]. 

The model (39) in the limit \ii = reduces to 

d 2 <p=- 7 ^-sm(qv 1 p 0V >). (42) 
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For later discussion we record here the mass of the soliton associated to this equation, 

8 

Correspondingly in the limit [13 = one has 



M mb = 7^-732 V^W- (43) 



d\ = -^sin(^A^) (44) 

with associated soliton mass 

M m = J^VfTi (45) 

Notice that other possibilities to perform the reduction of type (33) encounter 
some inconsistencies, e.g. the attempt to implement the reduction $1 = $3, $2 = 
q' <I>i implies 8\ 2 < which is a contradiction since are real numbers by defini- 
tion. The same inconsistency occurs when one tries to reduce the sl(3) GSG model 
to a three-frequency SG model. We expect that the three and higher frequency 
models [37] will be related to sl(N), N > 4, GSG models. 

In the following we will provide some kink solutions for particular set of para- 
meters. Consider 

v x = 1/2, Si = 6 2 = 1, v 2 = 1/2, 1/3 = 1/2 and q = 2, n = 2, m = 1 (46) 

which satisfy (38) and (40), respectively. This set of parameters provide the so-called 
double sine-Gordon model (DSG). Its potential — [4/xi(costt — l) + 2 / U3(cos<3? — 1)] has 
period 4ir and has extrema at <E> = 2-rrpi, and <!> = 47rp2±2cos~ 1 [l — |//i/(2/l*3)|] with 
Pi,p 2 £ the second extrema exists only if \^i/(2^)\ < 1. From the mathematical 
point of view the DSG model belongs to a class of theories with partial integrability 
[38]. Depending on the values of the parameters Pq, fii, ^3 the quantum field theory 
version of the DSG model presents a variety of physical effects, such as the decay 
of the false vacuum, a phase transition, confinement of the kinks and the resonance 
phenomenon due to unstable bound states of excited kink-antikink states (see [12] 
and references therein). The semi-classical spectrum of neutral particles in the DSG 
theory is investigated in [39]. Let us mention that the DSG model has recently been 
in the center of some controversy regarding the computation of its semiclassical 
spectrum, see [12, 13]. 

Interestingly the functions 1 

To = 1 + i dexp[y(x — vt)] + hexp[2j(x — vt)], 

T\ = 1 — i dexp[y(x — vt)] + hexp[2j(x — vt)], (47) 



1 These functions are obtained by adding the term exp[27(a; — vt)] to the relevant tau functions 
for one solitons used above. This procedure adds a new method of solving DSG which deserve 
further study. The multi-frequency SG equations can be solved through the Jacobi elliptic function 
expansion method, see e.g. [40]. 
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satisfy the equation (39) for the parameters (46) provided 
e-^' 2 = n/r 

7 2 = Mi + 2/X3, /l = _ Y' ei = ~ 1 
The general solution of this type can be written as 

"1 1 + h exp[27(x - vt)] 



$ := 4 arctan 



d exp[7(x — vt)] 



(48) 
(49) 



(50) 



2.5.1 DSG kink (h < 0, m > 0) 

For the choice of parameters h < 0, \ii > in (49) the equation (50) provides 



ip : = —arctan 
A) 



-2^V2 
d 



sinh[7^ (x — vt) + ao] 



IK 



±y/ Hi + 2^3, 
1. 



(51) 



a 



-\n\h\. 



This is the DSG 1-kink solution with mass 

Mi 



M K 



16 



2 1K 



1 + 



\/ 2 M3(Ml + 2 /U 3 ) 



VMi + 2/i 3 + \/2M3 v 
ln ( F= ) 



(52) 



Notice that in the limit iii — > the kink mass becomes M# = ^v/2u 3 , which 
is twice the soliton mass (43) of the model (42) for the parameters = 1/2, q = 2. 
Similarly, in the limit ^3 — > the kink mass becomes mSfp \[Wvi which is the soliton 
mass (45) of the model (44) for v\ = 1/2, q = 2; thus in this case the coupling 
constant is Po/2. As discussed above these solitons get in some sense "confined" 
inside the kink if the parameters satisfy \ii / 0. The 1-antikink is plotted in Fig. 
4. Moreover, the relevant baryon number associated to this DSG kink becomes 
B kmk = 4N c [6]. 

2.5.2 Bounce-like solution (h > 0, Hi < 0) 

For the parameters h > 0, Hi < one gets from (50) 



(p := —arctan 

Po 



2h 1 ' 2 
d 



cosh[7'(x — vt) + a' ] 



i = 2^3 - ImiI 



a 



1 



ln/i (53) 



This is the bounce-like solution and interpolates between the two vacuum values 
2ir and 4tt — 2arcos(l — |/xi/2/x 3 |) and then it comes back. Since 2ir is a false vacuum 
position this solution is not related to any stable particle in the quantum theory 
[12]. In Fig. 2 we plot this profile. 
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Figure 2: Bounce-like solution (Ifip) plotted for ^ 2 = -0.0000001, fi 3 = 
0.001, d = -2. 

3 Classical GSG as a reduced Toda model coupled to 
matter 

In this section we provide the algebraic construction of the sl(3) affine Toda model 
coupled to matter fields (ATM) and closely follows refs. [15, 34, 41] but the reduc- 
tion process to arrive at the classical GSG model is new. The previous treatments 
of the sl(3) ATM model used the symplectic and on-shell decoupling methods to 
unravel the classical GSG and generalized massive Thirring (GMT) dual theories 
describing the strong/weak coupling sectors of the ATM model [14, 15, 42]. The 
ATM model describes some scalars coupled to spinor (Dirac) fields in which the 
system of equations of motion has a local gauge symmetry. In this way one in- 
cludes the spinor sector in the discussion and conveniently gauge fixing the local 
symmetry by setting some spinor bilinears to constants we are able to decouple the 
scalar (Toda) fields from the spinors, the final result is a direct construction of the 
classical generalized sine-Gordon model (cGSG) involving only the scalar fields. In 
the spinor sector we are left with a system of equations in which the Dirac fields 
couple to the cGSG fields. 

The zero curvature condition (119) gives the following equations of motion [41] 

(54) 

im^-^^l + imle^-^fyl + im^e^^^l + m 2 e 3 '%55) 
my + ^^ R , -2d + ^ 2 L = mle^^l (56) 

(777 ^ 111 ^ \ 1 /2 
vm.<. J 



fP0a 

4i 

a = 1,2 
4~ 

-2d + ^ L = 
2d-^ R = 
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(7Th^~ 771^ \ 1/2 ■ 
-M) ^(VlVie^ 1 + ^e"*), (58) 

1 9 

(777 777 \ 1/2 
-*-±\ e \-$\tf R jfo +il) 2 L ip l R e i<tn ), (59) 



2<9_Vi = m^-^L 2cL^ = mje^^i, (60) 
-25 + & = m i e ^i^ + 2i(^^) 1/ %^(-^e- i *»-^e < *»), (61) 

1 ^ 

(777, , 777, , \ 1/2 . ~ 

) e"(ViV|e-^ + ^ e ^), (62) 

2<9_V% = mje" +i * 1 ^, -23+^1 = mje" - ** 3 ^, (63) 

-V 1 ) ^(VWIe^ - ^Vle^ 1 ), (64) 

2„ 



= 0, (65) 
where 0i = 26>i — 62, 4>2 = 2#2 — #i, 03 = #i + #2- Therefore, one has 

03 = 01 + 02 (66) 

The fields are considered to be in general complex fields. In order to define 
the classical generalized sine-Gordon model we will consider these fields to be real. 

Apart from the conformal invariance the above equations exhibit the (u (1)l) <8> 
(u (1)_r^ left-right local gauge symmetry 

a - ^ + e+(^+)+C(x-), a = 1,2 (67) 
? — > z7 ; 77 — > 77 (68) 

^ _> gi(l+75)S^(a; + )+j(l-75)S ! l(2;-) ^ (59) 
^ e -*(l+75)(S l + )( :E+ M(l-75)(H l _)(x-)^ ) ■ = 1)2)3; ( ? q) 

E\ = ±e± T 2&, S| = ±e± T 2&, El = ~l + E\. 

One can get global symmetries for £± = = constants. For a model defined by 
a Lagrangian these would imply the presence of two vector and two chiral conserved 
currents. However, it was found only half of such currents [34]. This is a consequence 
of the lack of a Lagrangian description for the sl(3)^ CATM in terms of the B and 

fields (see Appendix). So, the vector current 

3 

J^ = Y^ mffi-fiP (71) 
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and the chiral current 

3 

J^ = Y, "^'-/W + 2fy(mJ,0i + m^ 2 ) (72) 

are conserved 

d^ = 0, d li J 5 ' i = (73) 

The conformal symmetry is gauge fixed by setting 

r\ = const. (74) 

The off-critical model obtained in this way exhibits the vector and topological 
currents equivalence [41, 42] 

3 

J2 m^l^ 1 = e^d v {m\d x + m\d 2 ), m% = m\ + m%, m\ > 0. (75) 
j'=i 

Moreover, it has been shown that the soliton type solutions are in the orbit of 
the vacuum n = 0. 

In the next steps we implement the reduction process to get the cGSG model 
through a gauge fixing of the ATM theory. The local symmetries (67)- (70) can be 
gauge fixed through 

iipi'ip 3 = iAj = const.; ip 3 ry 5 i/j J = 0. (76) 
From the gauge fixing (76) one can write the following bilinears 

^Vl + ^ = 0, j = 1,2,3; (77) 

so, the eqs. (76) effectively comprises three gauge fixing conditions. 

It can be directly verified that the gauge fixing (76) preserves the currents con- 
servation laws (73), i.e. from the equations of motion (54)- (65) and the gauge fixing 
(76) together with (74) it is possible to obtain the currents conservation laws (73). 

Taking into account the constraints (76) in the scalar sector, eqs. (54), we arrive 
at the following system of equations (set r] = 0) 

d 2 6 1 = Mlsmfa + M^smfo, (78) 



/ , . Mill." .I/;'. .-,1 ",• = 

Define the fields ipi, tp 2 as 



d 2 6 2 = Mjsin^ + Mjsin^s, Mi > = 4A i mL i = 1,2,3. (79) 



Vl = a9 1 + b9 2 , a=-f d=-± 2 - 80 

3po^i^2 3/^0 ^1^2 

<P2 = c6 1 + d6 2 , b=-c= 2{ ll~ V2 \ vi,V2 el (81) 

3po^i^2 
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Then, the system of equations (78)-(79) written in terms of the fields 991,2 be- 
comes 

SPtp! = aM\ sin[/Vi(2^i - <p 2 )] + bM^ sm[(3 Q v 2 (2y 2 - <fi)] + 

(a + b)Ml sin/3 [(2^i - v 2 )^i + {2u 2 - ^2)], (82) 

d 2 tp 2 = cM\ sin[/3 ^i(2^i - ip 2 )] + dM\ sm\fov 2 {2ip 2 - <pi)] + 

(c + d)Ml sin/3 [(2^i - v 2 )y 1 + {2u 2 - u^)] (83) 

The system of equations above considered for real fields <pi 2 as well as for 
real parameters M^, a, b, c, d, (3q defines the classical generalized sine-Gordon model 
(cGSG). Notice that this classical version of the GSG model derived from the ATM 
theory is a submodel of the GSG model (10)-(11), defined in section 2, for the par- 
ticular parameter values r = 2vx ~" 2 , s = 2u2 ~ A Vl and the convenient identifications 
of the parameters in the coefficients of the sine functions of the both models. 

The following reduced models can be obtained from the system (82)-(83): 

i) SG submodels 

1.1) For v 2 = 2v\ one has = M 2 and the system 

w = 0, dV = M l^t sin A^m- 

1.2) For v\ = 2v 2 one has = and the system 
Vl = 0, d\ 2 = M^ sin/3 2^^2. 

1.3) For v 2 = v\ = v and 991 = ip 2 = 0a, (A = 1, 2), one gets the sub-models 
i.3a) M\ = M 2 ,M\ = 0, 3 2 (p x = aM\ sin fti/ft, 

i.3b) M\ = M 2 = 0, 3 2 <p 2 = aM] sm[3 o v0 2 . 

ii) DSG sub-model 

For v\ = v 2 and = M 2 one gets the sub- model ip± = ip 2 = ip, d 2 cp = 

aM^ sin (3 vi if + aM| sin 2/3 ^i f- 

The sub-models i.l)-i.2) each one contains the ordinary sine-Gordon model (SG) 
and they were considered in the subsections 2.1 and 2.2, respectively; the sub-model 
i.3) supports two SG models with different soliton masses which must correspond to 
the construction in subsection 2.3; and the ii) case defines the double sine-Gordon 
model (DSG) studied in subsection 2.5. Other meaningful reductions are possible 
arriving at either SG or DSG model. Notice that the reductions above are particular 
cases of the sub-models in subsections 2.1, 2.2, 2.3 and 2.5, respectively, for relevant 
parameter identifications. 

The spinor sector in view of the gauge fixing (76) can be parameterized conve- 
niently as 




(84) 



Therefore, in order to find the spinor field solutions one can solve the eqs. (56)- 
(64) for the fields Uj, v j for each solution given for the cGSG fields ip± i 2 of the system 
(82)-(83). 
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3.1 Physical solitons and kinks of the ATM model 

The main feature of the one 'solitons' constructed in [34] is that for each positive 
root of si (3) there corresponds one soliton species associated to the fields (f>±, 4>2, <f>3, 
respectively. The relevant solutions for the spinor fields together with the l-'solitons' 
satisfy the relationship (75). The class of 2-'soliton' solutions of sl(3) ATM obtained 
in [34] behave as follows: i) they are given by 6 species associated to the pair 
(aij, otj), i < j; i,j = 1, 2, 3; where the a's are the positive roots of si (3) Lie algebra. 
Each species (a^aj) solves the sl(2) ATM submodel 2 , ii) they satisfy the U(l) 
vector and topological currents equivalence (75). However, the possible kink type 
solutions associated in a non-local way to the spinor bilinears and the relevant gauge 
fixing of the local symmetry (67)-(70) have not been discussed in the literature. In 
order to consider the physical spectrum of solitons and study its properties, such as 
their masses and scattering time delays, it is mandatory to take into account these 
questions which are related to the counting of the true physical degrees of freedom 
of the theory. Therefore, one must consider the possible soliton type solutions 
associated to each spinor bilinear. The relation between this type of 'solitons', 
say 4>j, and their relevant fermion bilinears must be non-local as suggested by the 
equivalence equation (75). So, we may have soliton solutions of type 



At this stage one is able to enumerate the physical 1-soliton (1-antisoliton) spec- 
trum associated to the gauge fixed ATM model. In fact, we have three 'kinks' and 
their corresponding 'anti-kinks' associated to the fields <pi (i=l,2,3), and three kink 
and antikink pairs of type j = 1, 2, 3. Thus, we have six kink and their relevant 
antikink solutions, but in order to record the physical soliton and anti-soliton exci- 
tations one must take into account the four constraints (66) and (77). Therefore, 
we expect to find four pairs of soliton and anti-soliton physical excitations in the 
spectrum. This feature is nicely reproduced in the cGSG sector of the ATM model; 
in fact, in the last section we were able to write four usual sine-Gordon models as 
possible reductions of the cGSG model. Namely, one soliton associated to the fields 
ipi, (f2, respectively (subsections 2.1 and 2.2) and 1-solitons associated to the field 
fi = = fAt A = 1,2, respectively (subsection 2.3). In the 2-kink (2-antikink) 
sector a similar argument will provide us ten physical 2-solitons and their relevant 
2-antisoliton excitations, i.e. six pairs of 2-kink and 2-antikink solutions of type (f> 
and (f>, respectively, which give twenty four excitations, and taking into account the 
constraints (66) and (77) we are left with ten pairs of 2-solitons and 2-antisolitons. 
In fact, these ten 2-solitons correspond to the pairs we can form with the four species 

2 sl(2) ATM gauge unfixed 2— 'solitons' satisfy an analogous eq. to (75). Moreover, for ip real 
and ip — ±(V0* one nas i soliton-soliton SS, SS bounds and no SS (S =soliton, S=anti-soliton) 
bounds [43] associated to the field ip. 




(85) 
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of 1-solitons in all possible ways. The same argument holds for the corresponding 
ten 2-antisolitons. 

In this way the system (82)- (83) gives rise to a richer (anti)soliton spectrum 
and dynamics than the a field 'soliton' type solutions of the gauge unfixed model 
(54)-(64) found in [34]. Regarding this issue let us notice that in the procedure 
followed in ref. [34] the local symmetry (67)-(70) and the relevant gauge fixing has 
not been considered explicitly, therefore their 'solitons' do not correspond to the 
GSG solitons obtained above. 

Notice that the tau functions in section 2 possess the function j(x — vt) in their 
exponents, whereas the corresponding ones in the ATM theory have two times this 
function [34, 43]. This fact is reflected in the GSG soliton solutions which are two 
times the relevant solutions of the ATM model. It has been observed already in the 
sl(2) case that the 6 'soliton' of the gauge unfixed sl(2) ATM model (see eq. (2.22) 
of [43]) is half the soliton of the usual SG model. 



4 Topological charges, baryons as solitons and confine- 
ment 

In this section we will examine the vacuum configuration of the cGSG model and 
the equivalence between the U(l) spinor current and the topological current (75) 
in the gauge fixed model and verify that the charge associated to the U(l) current 
gets confined inside the solitons and kinks of the GSG model obtained in section 2. 

It is well known that in 1 + 1 dimensions the topological current is defined as 
j£ p ~ e^d^Q, where <3? is some scalar field. Therefore, the topological charge is 

(?top = J J® Q pdx ~ <£(+oo) — oo). In order to introduce a topological current 
we follow the construction adopted in Abelian affine Toda models, so we define the 
field 

o=l a 

where a a , a = 1,2, are the simple roots of sl(3). We then have that 6 a = (0|A O ), 
where A a are the fundamental weights of si (3) defined by the relation [44] 

2^f4 = S ab . (87) 
{a a \a a ) 

The fields <pj in the equations (54)- (64) written as the combinations (9\aj), j = 
1,2,3, where the a'^s are the positive roots of sl(3), are invariant under the trans- 
formation 



6^6 + 2ith or 0j-»0j + 27r(/i|ay), (88) 

{a a \a a ) 



a* - E^rrH (89) 



z 
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where \i is a weight vector of sl(S), these vectors satisfy (fJ>\ctj) G 2 and form an 
infinite discrete lattice called the weight lattice [44]. However, this weight lattice 
does not constitute the vacuum configurations of the ATM model , since in the 
model described by (54)- (65) for any constants 9^ and r/°) 

^ = ^ = 0, 9 a = df\ 7] = r/°\ v = -mV^x+aT (90) 

is a vacuum configuration. 

We will see that the topological charges of the physical one-soliton solutions of 
(54)-(65) which are associated to the new fields <p a , a = 1,2, of the cGSG model 
(82)-(83) lie on a modified lattice which is related to the weight lattice by re-scaling 
the weight vectors. In fact, the eqs. of motion (82)-(83) for the field defined by 
if = Yla=i ^r^a, such that ip a = (<p|A a ), are invariant under the transformation 

2 

— i — \>*» gZ - (91) 
Po ^-f v a (a a \a a ) 

a=L 

So, the vacuum configuration is formed by an infinite discrete lattice related 
to the usual weight lattice by the relevant re-scaling of the fundamental weights 
A a — > ^A a . The vacuum lattice can be given by the points in the plane ipi x ip 2 

2vr 2 qi q 2 q-y 2q 2 

(<Pi , <P2) = + — ,— + ), 9«6 2. (92) 

3po v\ v 2 v\ v 2 

In fact, this lattice is related to one in eq. (9) through appropriate parameter 
identifications. We shall define the topological current and charge, respectively, as 

J top = l^^' Q top = / dxJ^p = §k(+oo) - tfoo)]. (93) 

Taking into account the cGSG fields (82)-(83) and the spinor parameterizations 
(84) the currents equivalence (75) of the ATM model takes the form 

3 

m^W = e^d^ y,! + Cl <P2), (94) 

j=i 

where = 0QViv 2 {m^d + m^b), (| = &QV\v 2 {m^a - m^b) and the spinors are 
understood to be written in terms of the fields Uj and Vj of (84). 

Notice that the topological current in (94) is the projection of (93) onto the 

vector f (^i + C|A 2 ). 

As mentioned in section 3 the gauge fixing (76) preserves the currents conser- 
vation laws (73). Moreover, the cGSG model was defined for the off critical ATM 
model obtained after setting r\ = const. = 0. So, for the gauge fixed model it is 
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expected to hold the currents equivalence relation (75) written for the spinor pa- 
rameterizations Uj,Vj and the fields <pi t 2 as is presented in eq. (94). Therefore, in 
order to verify the U(l) current confinement it is not necessary to find the explicit 
solutions for the spinor fields. In fact, one has that the current components are 
given by relevant partial derivatives of the linear combinations of the field solutions, 
i.e. J° = E -=i m^ 7 V = d x (Cl ¥>i + and J 1 = £* =1 m^W = 

—dt{Clp ¥>i + C$ V^)- In particular the current components J , J 1 and their associated 
scalar field solutions are depicted in Figs. 3 and 4, respectively, for antisoliton and 
antikink solutions. 

It is clear that the charge density related to this U(l) current can only take 
significant values on those regions where the x— derivative of the fields </?i ; 2 are non- 
vanishing. That is one expects to happen with the bag model like confinement 
mechanism in quantum chromodynamics (QCD). As we have seen the soliton and 
kink solutions of the GSG theory are localized in space, in the sense that the scalar 
fields interpolate between the relevant vacua in a limited region of space with a size 
determined by the soliton masses. The spinor U(l) current gets the contributions 
from all the three spinor flavors. Moreover, from the equations of motion (56)-(64) 
one can obtain nontrivial spinor solutions different from vacuum (90) for each set 
of scalar field solutions (fi,ip2- For example, the solution ipi =soliton, 1^2 = in 
section 2.1 implies fa = (pi, fa = —'-Pi-, fa = Q which substituting into the spinor 
equations of motion (56)- (64) will give nontrivial spinor field solutions. Therefore, 
the ATM model of section 3 can be considered as a multiflavor generalization of the 
two-dimensional hadron model proposed in [30, 31]. In the last reference a scalar 
field is coupled to a spinor such that the DSG kink arises as a model for hadron and 
the quark field is confined inside the bag. 
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Figure 3: 1-antisoliton and confined current J M . The solid curve is the 1- 
antisoliton ip), the dashdotted curve is J° and the curve with losangles is 
J 1 . For t = 1, /ii = n% = 1, d = 1.5, v = 0.05, /3o = 0.5, m\p = = l,ui = 
1, tfi = 1, 6 2 = 2. 

5 Qualitons or quark solitons in two-dimensional QCD 

Several properties of the ATM model deserve careful consideration in view of the 
relationships with two-dimensional QCD. In particular, it has been shown that the 
sl{2) ATM model describes the low-energy spectrum of QCD2 (1 flavor and N c col- 
ors) [19]. In the context of bosonized QCD2 the appearance of soliton solutions that 
have the quantum numbers of quarks as constituent of hadrons has been considered 
[8]. So, one can inquire about these type of quark solitons in the context of the 
ATM model description of QCD2. Since the ATM model describes the low-energy 
effective action in the strong coupling limit of QCQ2, in order to disentangle the 
quark solitons one needs to restore, in some way, the heavy fields, i.e. the fields 
associated to the color degrees of freedom. For simplicity we choose the s/(2) case 
in the following developments. 

The Lagrangian of the si (2) ATM model is defined by [42, 43, 45] 

= ^ipd^ip + i^d^-m^p^e 2 ^^^, (95) 

where k = (k G Z), ip is a real field, is a mass parameter, and ip is a 
Dirac spinor. Notice that rp = -0 T 7o- We shall take ip = e^ijj* [43], where is a 
real dimensionless constant. The conformal version (CATM) of (95) has been con- 
structed in [41]. The integrability properties and the reduction processes: WZNW^ 
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Figure 4: DSG kink solution and confined current J^. The curve with 
losangles is the antikink (fyf), the dashdotted line is J°, the solid curve is 
J 1 . For t = l,po = 10 8 , rnf = m = -0.0000001, /x 3 := 0.001, d = 2,Sj.= 
6 2 = l,vi = 1/2. 

CATM -> ATM -> sine-Gordon(SG) + free field, have been considered [42, 43, 45] 
. The sl(n) ATM exhibits a generalized sine-Gordon/massive Thirring correspon- 
dence [14]. Moreover, (95) exhibits mass generation despite chiral symmetry [46] 
and confinement of fermions in a self-generated potential [43, 30]. 

The Lagrangian is invariant under ip — > (p + mr, thus the topological charge, 
Qtonol = j dxj°, = ^e^d^ip, can assume nontrivial values. A reduction is 
performed imposing the constraint 

-L e /^ <p = I^ 7 />, (96) 

Ztt tt 

where = ip^^ip is the U(l) Noether current. In fact, the soliton type solutions 
satisfy this relationship [43]. 

The Eq. (96) implies ip'tp ~ d x ip, thus the Dirac field is confined to live in 
regions where the field ip is not constant. The 1(2)— soliton(s) solution(s) for cp and 
ip are of the sine-Gordon (SG) and massive Thirring (MT) types, respectively; they 
satisfy (96) for |e^| = 1, and so are solutions of the reduced model [43]. Similar 
results hold in sl(n) ATM [34, 14]. 

The equivalence (96) for multisolitons describes, ip = ipjy [Qtopol = -^ sl S n ( e '</')] 
and N— solitons of the SG and MT type, respectively. Asymptotically one can 
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write 

1 N 1 

Z7T z — ' 7T 

a=l 

where the ip a 7 s are the solutions for the individual localized lowest energy fermion 
states. In fact, (97) encodes the classical SG/MT correspondence [47]. Thus, the 
ATM model can accomodate iV c = N— fermion confined states with internal 'color' 
index a [30]. 

In order to gain insight into the QCD2 origin of the ip a fields let us write the 
'mass term' in the multifermion sector of ATM theory as [19] 

i, a e 2 ^5 r = e 2i<p + ^ Lae -2up m (98) 

The ATM mass term in the multifermion sector, Eq. (98), must be compared to the 
corresponding term in the bosonized QCD2 in order to identify the fields related to 
the flavor and color degrees of freedom, respectively. 

Therefore the total chiral invariant Lagragian including the kinetic terms for the 
quark fields becomes 

±£ = \d^d^ + ie^J2(^d^ a -m^^ a e 2i ^^ a ). (99) 

a 

Although the QCD color degrees of freedom have a non-abelian symmetry we 
use abelian bosonization techniques in order to bosonize the fermions. This will be 
sufficient in order to reproduce various properties of the effective QCD2 Lagrangian 
in this regime as presented in ref. [8]. So, let us introduce new boson field represen- 
tations of the fermion bilinears as [48] 



i:^ a rd^ a : = - — (d u 4> a ) 2 (100) 
:^ a (l±75^ a : = _^ :e (±-<M : , :^ a rf^ a :=--^d v tf> a , (101) 

7T 7T 

where c = \ exp (7), /x is an infrared regulator and a a real parameter. 

In order to compare to the related QCD Lagrangian describing the regime m q >> 
e c [8], which does not possess an exact chiral symmetry, we must introduce some 
chiral symmetry breaking terms in the Lagrangian (99). The most direct program 
for accomplishing this is simply to include certain chiral breaking terms in the 
bosonized version of the ATM+color model given in (99) in the form of 

C bos = fo„pfflV + + fcmW cos(2y + 2o^a) - m»#} - 

a 

"W 2 - ^2 m ab(f>a4>b ~ ^2 m 0aV4>a (102) 
a<b a 
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Notice that we have included certain bilinear terms in the scalar fields as the sym- 
metry breaking terms. Define the fields \a and as 

Xa=-^a + ^ * = ^-_g Xa ), d=- + ^^. (103) 

So, providing the relationships 

m a = const, rrigja = m^ a = const, (a < b), niQ a = const, Va, (104) 

mab ~ 2m N c ' Wa " 4m + iV c ' dc " p ' (1U5) 

^ = (106) 
the Lagrangian (102) becomes 

C bos = ^(^) 2 -^(l + ^)<I> 2 + j:{^W + 2M 2 (cos/3 Xa )}- 

a 

^ E - 2e c(^) E X« + 4 4 E *»X> (^7) 



a<b a c a <b 



where 



r = *. = "* = WF. - 1) <108) 

The model (107) except the a < b kinetic (the first term of the second line in 
(107)) and the $ terms reproduces the QCD2 bosonized Lagrangian (in the regime 
m q » e c ) presented in [8]. Notice that the $ field completely decouples from the 
rest of the fields. Moreover, in the opposite limit, i.e. the strong coupling regime 
and large N limit we can verify that this field becomes a free massless field [19]. 

Besides, the low-energy spectrum of QCD2 has ben studied by means of abelian 
[49] and non-abelian bosonizations [50, 2]. In this limit the baryons of QCD2 are 
sine-Gordon solitons [2]. In the large N limit approach (weak e and small m q ) the 
SG theory also emerges [51]. 

The question of confinement of the "color" degrees of freedom associated to the 
field V' in the ATM model by computing the string tension has been presented in 
[19]. In the fundamental representation of the quarks it has been taken = m q 
and k = 2N/-K. Then from (108) one has \e f \ = 2y/ £_ N - 

Following [8] we define the baryon number as 

, N ' 

B 



^E [xfc(+°o)-Xfc(-°o)] ( 109 ) 



7T 

k=l 
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We seek for solutions of the field equations of motion in the static case 

N c N c 
xl-iM^senV^Xa + P^x'l-lel^Xa-j^Y.Xb) = °> ( U0 ) 

b>a C 6=1 

p = [N - 1 + 4iry/N{N - I)] -1 . (Ill) 

Depending on the boundary conditions for the fields Xfc(±°o) we may have cer- 
tain nucleon states with B = kN c , k £ 7L. (the baryon number is normalized to be 
N c for the nucleon) or some quark solitons (B = n, n =integer non-multiple of N c ). 
These type of solutions can be discussed by analyzing the field equations for the 
static case [8]. In the low energy and strong coupling limit (e c >> m q ) the nucleon 
states (baryons and multibaryon) are described by the generalized sine-Gordon soli- 
tons (see [6] and references therein), as can be inferred from the form of the eqs. of 
motion (110) in this limit. Whereas, the quark solitons exist for a sufficiently heavy 
quark m q , but have infinite energy, corresponding to a string carrying the non-singlet 
color flux off to spatial infinity, i.e. they exist in the opposite limit m q » e c . These 
quark soliton solutions disappear when the meson mass parameter M is reduced to 
become comparable to the gauge coupling strength e c (it has the dimension of mass 
in QCD2). Let us search for solutions such that Xa(— oo) = for all a. Then, at 
x = +oo one has 




(112) 



The eq. (112) becomes the same as the one presented in [8] describing the boundary 
condition at x = +oo. If we assume Xa(oo) = x f° r an a > one nas that x(°°) = \\f^ n ^ 
and B = ^nN c . But, in order to have positive eigenvalues of the squared mass matrix 

9x dxb ' We mus ^ nave even n, and thus integer baryon number B = kN c (baryons 
and multibaryons) . 

Following [8] , in the search for quark solitons let us first concentrate on the case 
N c = 2. Then, eq. (112) can be written as 

sinA/47rxi(°o) = -eV7r[xi(°°) -X2(oo)], (113) 
sinv / 4vrX2(oo) = -ev^^tx)) - Xi(oo)], (114) 

e 2 

where e = 2n M^ ■ We may have non-baryonic solitons with B = n for odd values of 
n (the quarks correspond to n = 1). For e << 1 we can find a series of solutions 
with positive second derivative matrix. This solution satisfies 

X2(oo) = -X2(oo) + n\pK, (115) 

which together with (113)- (114) provides 

sin\/4vrxi = — e(V47TXi — n7r )- (H6) 
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Figure 5: The axis is defined as x = \/47r\;i (00) . Comparison of the left-and 
right-hand sides of the soliton eq. (116), corresponding to a quark soliton 
with n = 1, for e = 0.08 (dotted line) and e = 0.4 (solid line). There are no 
solutions for e > 1 (dot-dashed line, drawn for e = 1.5). 



Let us define £ = v4jt[x(oo) — ^n^/ir], then one has that the solutions are 

f (tt- e)(2i) for n even 
Q \ (7r-e)(2i + l) for n odd ' 1 j 

in the limit where (le << j). The solutions (117) correspond to excitations of 
"colored" states and have infinite energy, with classical string tension 



7r (2Z) 2 for n even 
vr e 2 (2Z + l) 2 for n odd 



2 /or 1 1 ^2 ^ „ j j 

(118) 



The single constituent quark soliton corresponds to n = 1, (2/ + 1) = 1. Thus, we 
have shown that QCD2 has quark soliton solutions if the quark mass is sufficiently 
large. These quark solitons disappear when the quark mass m q is reduced until the 
meson mass M becomes comparable to the dimensional gauge coupling strength e c . 
The above picture can be directly generalized for any N c , see more details in [8]. 



6 Discussion 



The generalized sine-Gordon model GSG (10)-(11) provides a variety of solitons, 
kinks and bounce type solutions. The appearance of the non-integrable double sine- 
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Gordon model as a sub-model of the GSG model suggests that this model is a non- 
integrable theory for the arbitrary set of values of the parameter space. However, a 
subset of values in parameter space determine some reduced sub-models which are 
integrable, e.g. the sine-Gordon submodels of subsections 2.1, 2.2 and 2.3. 

In connection to the ATM spinors it was suggested that they are confined inside 
the GSG solitons and kinks since the gauge fixing procedure does not alter the 
U(l) and topological currents equivalence (75). Then, in order to observe the bag 
model confinement mechanism it is not necessary to solve for the spinor fields since 
it naturally arises from the currents equivalence relation. In this way our model 
presents a bag model like confinement mechanism as is expected in QCD. 

The (generalized) massive Thirring model (GMT) is bosonized to the GSG model 
[17], therefore, in view of the solitons and kinks found above as solutions of the GSG 
model we expect that the spectrum of the GMT model will contain 4 solitons and 
their relevant anti-solitons, as well as the kink and antikink excitations. The GMT 
Lagrangian describes three flavor massive spinors with current-current interactions 
among themselves. So, the total number of solitons which appear in the bosonized 
sector suggests that the additional soliton (fermion) is formed due to the interactions 
between the currents in the GMT sector. However, in subsection 2.3 the soliton 
masses M3 and M4 become the same for the case n\ = [12, consequently, for this 
case we have just three solitons in the GSG spectrum, i.e., the ones with masses 
Mi, M2 (subsections 2.1-2.2 ) and M3 = M4 (subsection 2.3), which will correspond 
in this case to each fermion flavor of the GMT model. Moreover, the s/(3) GSG 
model potential (6) has the same structure as the effective Lagrangian of the massive 
Schwinger model with Nf = 3 fermions, for a convenient value of the vacuum angle 
9. The multiflavor Schwinger model resembles with four-dimensional QCD in many 
respects (see e.g. [52] and references therein). 

The sl{n) ATM models may be relevant in the construction of the low-energy ef- 
fective theories of multiflavor QCD2 with the dynamical fermions in the fundamental 
and adjoint representations. Notice that in these models the Noether and topologi- 
cal currents and the generalized sine-Gordon/massive Thirring models equivalences 
take place at the classical [15, 42] and quantum mechanical level [17, 43]. 

The interest in baryons with exotic quantum numbers has recently been sti- 
mulated by various reports of baryons composed by four quarks and an antiquark. 
The existence of these baryons cannot yet be regarded as confirmed, however, re- 
ports of their existence have stimulated new investigations about baryon structure 
(see e.g. [53] and references therein). Recently, the spectrum of exotic baryons in 
QCD2, with SU(Nf) flavor symmetry, has been discussed providing strong support 
to the chiral-soliton picture for the structure of normal and exotic baryons in four 
dimensions [6, 54]. The new puzzles in non-perturbative QCD are related to sys- 
tems with unequal quark masses, so the QCD2 calculation must take into account 
the SU (N f)-breakmg mass effects, i.e. for Nf = 3 it must be m s / m u ^- So, in 
view of our results above, the properties of the GSG and the ATM theories may 
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find some applications in the study of mass splitting of baryons in QCD2 and the 
understanding of the internal structure of baryons. Regarding this line of research, 
it has been shown that the GSG model describes the low-energy spectrum of normal 
and exotic baryons in QCD2 with unequal quark mass parameters [6]. 

Finally, we have considered the quark soliton (qualiton) solutions of QCD2 in 
the regime e c « m q . In this context the role played by the sl(2) ATM model is 
clarified. In fact, the qualitons arise if the color degrees of freedom are restored by 
coupling them to the Toda field and convenient boundary conditions are imposed 
on the fields. So, we have shown that the sl(2) ATM model becomes a low-energy 
effective lagrangian describing the quark confinement mechanism in QCD2. The 
equivalence between the Noether and topological currents (96) is a crucial property 
of the ATM model in order to provide the confinement mechanism. This picture 
can be directly generalized to any number of flavors Nf since a relationship analog 
to (96) holds in that case, e.g. the Nf = 3 case is presented in (75). 
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A The zero-curvature formulation of the ATM model 

We summarize the zero-curvature formulation of the si (3) ATM model [14, 15, 34]. 
Consider the zero curvature condition 

d+A_ - d-A + + [A+, A-] = 0. (119) 

The potentials take the form 

A + = -BF+B' 1 , A- = -d-BB~ l + F~, (120) 

with 

F+ = F+ + F+, F~ = Ff + F 2 ", (121) 
where B and F^ 1 contain the fields of the model 

Ft = V ^^< + V ^V'|< + V^4V ; ^-a3 ) (!22) 

Ff = y^^< + V^^-ai + V^4^-^' ( 123 ) 

F{ = ^^lE-l + ^^\E\ l + ^^lE Q _ a2 , (124) 
F 2 = ^^\F~l + ^^\E~l + y^4V%£° a3 , (125) 

£> = e i9iH ( l+id 2 H% e »C e r)Q ppal = jjgVC e r)Q ppal ^ (126) 
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E™.,H™, H 2 and C (i = 1, 2, 3; n = 0, ±1) are some generators of sl(3)^; Q ppa i 
being the principal gradation operator. The commutation relations for an affine Lie 
algebra in the Chevalley basis are 

[B?, m = mC^K ab 8 m+n ,o (127) 
[K,El a ]=±K aa E^+ n (128) 

K, E -a] = £ l a aK l+n + ^mCS m+nfi (129) 

a=l 

[E%, E%] = e(a, (3)E™+£; if a + [5 is a root (130) 
[D,E n a ] = nE n a , [D,H n a ]=nH n a . (131) 

where -fT aa = 2a.a a /a^ = n^K ba , with ra" and being the integers in the expansions 
a = n"a a and a/a 2 = /"a a /a 2 , and e(a,{3) the relevant structure constants. 

Take K\\ = K 22 = 2 and i^i2 = K 2 i = — 1 as the Cartan matrix elements 
of the simple Lie algebra sl(3). Denoting by ol\ and a 2 the simple roots and the 
highest one by ip(= u\ + a 2 ), one has if = l(a = 1, 2), and K^i = K^ 2 = 1- Take 
e(a,f3) = -e(-a,-f3), e lj2 = e{a 1 ,a 2 ) = 1, e_i j3 = e(-ai,V>) = 1 and e_ 2 ,3 = 
e(-a 2 ,i/j) = -1. 

One has Q ppa / = ^a=i s aKi-H + 3-D, where are the fundamental co-weights 
of sl(3), and the principal gradation vector is s = (1, 1, 1) [55]. 
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